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Kagome lattice is a two-dimensional network of corner-sharing triangles and is often associated 
with geometrical frustration. In particular, the frustrated coupling between waveguide modes in a 
kagome array leads to a dispersionless flat band consisting of spatially localized modes. Here we 
propose a complex photonic lattice by placing T^T-symmetric dimers at the kagome lattice points. 
Each dimer corresponds to a pair of strongly coupled waveguides. With balanced arrangement of 
gain and loss on individual dimers, the system exhibits a PT-symmetric phase for finite gain/loss 
parameter up to a critical value. The beam evolution in this complex kagome waveguide array 
exhibits a novel oscillatory rotation of optical power along the propagation distance. Long-lived 
local chiral structures originating from the nearly flat bands of the kagome structure are observed 
when the lattice is subject to a narrow beam excitation. 


Photonic lattices composed of balanced gain and loss 
waveguides urn have attracted considerable attention 
due of their potential applications in optical beam en¬ 
gineering and image processing mzi. These photonic 
lattices belong to a larger class of intriguing metamate¬ 
rials that exhibit the parity-time {VT) symmetry [SUTO] . 
The notion of T^T-symmetry was originally introduced 
as a complex extension of quantum mechanics m. It 
was shown that non-Hermitian Hamiltonians commut¬ 
ing with the combined VT operator may exhibit entirely 
real energy spectra. In optics, this combined symme¬ 
try requires that the complex refractive index obeys the 
condition n(r) = n*(—r). The system is then in the so- 
called exact phase with real propagation constants for 
a small gain and loss coefficient. As this coefficient in¬ 
creases, the system undergoes a transition into a phase 
with broken T^T-symmetry. Although optical realiza¬ 
tion of T^T-symmetric Hamiltonians is usually based on 
the formal analogies between Schrodingier and paraxial 
Helmholtz equations, it was demonstrated recently that 
T^T-symmetric optical phases exist in regimes well be¬ 
yond the paraxial approximation [12]. 

Motivated by recent experimental developments [U 
[3], T^T-symmetric photonic lattices have been theoret¬ 
ically studied in various geometries, including the one¬ 
dimensional (ID) chain jH [13] and ladder [14], as well 
as 2D square [4], honeycomb na ig, and triangular 
lattices m- Several intriguing features have been ob¬ 
served including, among others, power oscillation, dou¬ 
ble diffraction, and a new type of conical diffraction. In¬ 
clusion of nonlinearity also leads to several new soliton 
solutions [ITHlQ]. Another lattice structure of both fun¬ 
damental interest and practical application is the kagome 
lattice. The interest in it is partly due to the optical flat 
bands that have been recently demonstrated in the so- 
called Lieb lattice mi ng. In both kagome and Lieb 
lattices, strong geometrical frustration leads to a degen¬ 
erate flat band in their coupled-mode spectra. In par¬ 


ticular, kagome is one of the most extensively studied 
structure in highly frustrated magnets [23]. The eigen- 
modes of the flat band are highly localized spatially even 
in the linear regime. Diffraction-free transmission can be 
achieved in such lattices by decomposing a given image 
into the local modes [24]. The significantly reduced dis¬ 
persion of the flat band also indicates enhanced nonlinear 
effects [25[[26]. In addition to the flat band, the kagome 
spectrum also contains several topologically nontrivial 
band-crossing points; gapping out these points could lead 
to novel optical topological phases [27] . 

Our goal in this paper is to present a T^T-symmetric 
kagome lattice which can be used to explore the intrigu¬ 
ing interplay of flat band and novel phenomena due to 
complex potentials. However, a straightforward arrange¬ 
ment of gain and loss in a kagome lattice fails to pro¬ 
duce a T^T-symmetric phase. Our strategy is to build 
the complex waveguide array based on T^T-symmetric 
dimers, which have been shown to play an important 
role in realizing an exact T^T-symmetric phase in other 
structures [THITQ]. A T^T-symmetric dimer is a pair of 
strongly coupled waveguides with balanced gain and loss 
Tiq. In fact, optical lattices built from such dimers pos¬ 
sess a local VT symmetry associated with each dimer 
and are rather robust against disorder M- On the other 
hand, the dimerization also significantly lowers the lattice 
symmetry to no more than a mirror reflection. Indeed, 
in order to realize the exact VT phase, the Cs symmetry 
is lost in the complex triangular and honeycomb lattices 
proposed in Refs. [IMI]. 

In this paper, we propose and study a kagome photonic 
lattice built from T^T-symmetric dimers. The lattice 
structure can be viewed as two interpenetrating twisted 
kagome lattices and retains a full Cs point group symme¬ 
try. We show that an exact T^T-symmetric phase exists 
for gain/low parameter smaller than a critical value 7c. 
More importantly, the band structure with 7 < 7c con¬ 
tains two nearly flat bands inherited from the kagome 
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FIG. 1: (Color online) kagome photonic lattice consisting of V'T- 
symmetric dimers. The blue and red circles denote waveguides 
with a gain and a loss parameter izy, respectively, n is the intra¬ 
dimer coupling, while t and t' denote the nearest-neighbor (NN) 
couplings between same and different sublattices, respectively. 


tight-binding spectrum. The T^T-symmetry breaking 
transition occurs when the two nearly flat bands collapse 
at the Brillouin zone center. Our linear beam dynam¬ 
ics simulations uncover interesting oscillatory rotations 
of optical power along the propagation distance. In par¬ 
ticular, we find long-lived local chiral structures consist¬ 
ing of quasi-local modes when the lattice is subject to 
a narrow beam excitation. The non-reciprocal nature of 
the chiral beam transmission is also discussed. 

The complex kagome lattice depicted in Fig. con¬ 
sists of dimers placed at vertices of each triangle. Since 
the kagome-lattice sites correspond to bonds in the dual 
honeycomb lattice, the same structure can be obtained 
by placing dimers on the bonds of the honeycomb lat¬ 
tice, with dimer orientations perpendicular to the bonds. 
Each dimer consists of two types of waveguides: type 
A is made of gain material, whereas type B exhibits an 
equal amount of loss. Interestingly, the A sublattice itself 
forms a twisted kagome lattice, while the B waveguides 
are arranged in another twisted kagome with opposite 
chirality. We assume that each waveguide supports only 
one mode, while light is transferred between neighbor¬ 
ing waveguides through optical tunneling. In the tight- 
binding description [28l, the diffraction dynamics of the 
optical field amplitude = (a^, 6 ^)^ at the n-th dimer 
evolves according to the following coupled-mode equa¬ 
tion: 

da 

m 

db 

i-^ = -i'f bn + Kan + '^[tbjn + {f ± S) am] ■ (1) 

m 

Here 7 is the gain/loss parameter, n is the dominant 
intra-dimer coupling, t is the nearest-neighbor (NN) cou¬ 
pling constant between waveguides of same sublattices. 
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FIG. 2: (Color online) Band structure of the kagome dimer lattice 
along high-symmetry directions for varying gain and loss parameter 
7 . The coupling constants (in units of the NN coupling t) are k = 3, 
t' = 1.1, and (5 = 0.1. The critical 7 c = 0.8. (a) and (b) show the 
(completely) real spectra when 7 = 0.7 and 7 = 7c, respectively, 
(c) and (d) show the real and imaginary parts, respectively, of the 
eigen-spectrum for 7 = 1 . 2 . 


t' ± J denote the two distinct NN couplings between dif¬ 
ferent sublattices; see Fig. The summation above is 
restricted to the NN pairs. 

We next examine the dispersion relations of the com¬ 
plex kagome lattice. To this end, we note that the kagome 
site index can be represented as n = ( 5 ,r^), where 
s = 1, 2, 3 indicates the three sublattices of kagome, and 
Tn denotes the Bravais lattice point of the underlying 
triangular lattice. After introducing a Fourier transform 
for field amplitude = (I/VN) exp(ik-rn) and 

a similar expression for 6 ^, the coupled-mode equations 
become id'^i^jdz = Hk ^k, where 


Hk 


A\^ Sk \ 

^k ‘^k / 


and the two 3x3 matrices A and B are 


( 2 ) 


( iy tcs tc2 \ 

tcs ij tci , (3) 

tc2 tci iy ) 


( K t'c^ — iSss t'c2 — iSs2 \ 
t'c^ + i6ss K t'ci -h iSsi I . (4) 

t'c2 + iSs2 t'ci — idsi K J 

For convenience, we have introduced factors q = 2 cos(k' 
di) and Si = 2sin(k • d^) (i = 1,2,3), the three vectors 
di ,2 = ^^( 1 , Ta/ 3)/4, ds = a(l/ 2 , 0 ) connect the neigh¬ 
boring sites on kagome, and a is the lattice constant. 
The dispersion relations of the eigenmodes are given by 
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FIG. 3: (Color online) Snapshots of the intensity profile |^nP 
under wide beam excitations for (a) a passive array (7 = 0 ) and 
(b) a complex array with 7 = 0.75; the width of the incident beam 
is wq = 5.5a,. The lower panels show the radial distribution of 
optical power for (c) a passive and (d) a complex 7 = 0.75 kagome 
array. 

the stationary solutions ^ ex.p{if3z), where the prop¬ 
agation constant f3 = /3(k) is the eigenvalue of 'Hk- 

As is well established in other non-Hermitian systems, 
the VT symmetry is a necessary but not a sufficient con¬ 
dition for the reality of the eigenvalues By using spec¬ 
tral techniques we numerically determined the existence 
of a VT threshold 7 c, below which the propagation con¬ 
stants of all bands and wave vectors are real. The critical 
7 c depends on coupling constants k and t', to be discussed 
below. Above this threshold, the system undergoes a 
transition into a phase with partially complex eigenval¬ 
ues /3. Fig.j^a) shows the band structure of the kagome 
lattice when 7 < 7c- This spectrum inherits several fea¬ 
tures characteristic of the kagome tight-binding model; it 
can be viewed as two (scaled) copies of the kagome spec¬ 
trum. First, there are two pairs of Dirac points located at 
the corners of the Brillouin zone [25]. Second, there are 
two quadratic band crossing points at the zone center. 
These quadratic crossing points are topologically non¬ 
trivial in the sense that each of them carries a 27r Berry 
flux [27j. Finally, there are two nearly flat bands (the 3rd 
and 4th) due to geometrical frustration; the weak disper¬ 
sion is caused by the coupling between the two twisted 
kagome sublattices A and B. 

The two nearly flat bands are separated by a finite gap 
A as shown in Fig. [^a). With increasing gain/loss pa¬ 
rameter, the closing of the gap at k = 0 coincides with 
the T^T-symmetry breaking transition. Analytical calcu¬ 
lation finds a gap Ap = 2 ^/{n — 2t'Y ~ 7 ^ a,t the F point. 
Setting Ap = 0 gives the critical value for the VT tran- 



FIG. 4: (Color online) The oscillatory power rotations in a W- 
symmetric kagome lattice (7 = 0.78) for (a) wide and (b) narrow 
beam excitations. 


sition: 7 c = —2t', which is independent of t and 5. The 

two topological quadratic band-crossing points merge at 
the Brillouin zone center when 7 = 7c; see Fig. ib). 
The VT phase transition here exhibits all the character¬ 
istics of an exceptional-point singularity. Other than the 
coalescing of eigenmodes, signaling a collapsed Hilbert 
space, we have also confirmed numerically a divergent 
Petermann factor mED] when 7 ^ 7c. 

The complex guiding potentials also have profound ef¬ 
fects on linear beam dynamics of the kagome waveguide 
arrays. By exciting the T^T-symmetric lattice with a 
Gaussian beam at 2 ; = 0, we numerically integrate Eq. Q 
for a large system with up to 6 x 10^ waveguides. We 
first consider the situation when the array is excited by 
a wide beam of width wq = 5.5a; the center of the beam 
coincides with one of the triangle center. For a passive 
kagome lattice, the beam evolution follows the standard 
diffraction pattern as shown in Fig. [^a). The envelop 
of the optical field remains circularly symmetric with re¬ 
spect to the beam center, while its amplitude slowly de¬ 
cays with distance; see Fig. [^c). On the other hand, 
the intensity profile. Fig. |^b), of a complex array with 
7 = 0.75 exhibits a rather inhomogeneous power distribu¬ 
tion within individual dimers. As the optical power oscil¬ 
lates in the kagome lattice, the radial profile of |Tp shows 
multiple branches as shown in Fig.|^d). Detailed exam¬ 
ination reveals that the optical field rotates periodically 
about the beam center. To quantify this oscillatory rota¬ 
tion, we first compute the angular Fourier components of 
the intensity profile: Tm = exp(im 6 >n )5 where 

m is an integer, and 0^ = arctan(^^/x^) is the angular 
coordinate of the nth waveguide. The fact that the opti¬ 
cal field maintains a Cs symmetry indicates that the first 
nontrivial higher harmonic is m = 3. We can then define 
an angle ©3 = arctan[Im(F 3 )/Re(F 3 )] to measure the 
shift of the optical power from the Cs symmetric axis of 
the lattice. For a passive lattice, this angle remains zero 
throughout the beam propagation. On the contrary, ©3 
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FIG. 5: (Color online) Snapshots of the quasi-stationary local 
structure under narrow beam excitations (wq = a) for (a), (b) 
a passive array (7 = 0) and (c), (d) a complex array with 7 = 0.75. 


oscillates with z in the exact T^T-symmetric phase of the 
kagome lattice; see Fig. |^a). This result is reminiscent 
of the power oscillation phenomena of a single VT dimer. 

An input beam with a narrow width excites a broad 
range of eigenmodes including those in the nearly flat 
bands. After the optical power spreads over the lattice 
through the dispersive modes, a long-lived local structure 
remains in the central region of the input beam. These 
local structures are quasi-stationary. For example, in a 
passive array, the intensity profile oscillates between the 
two patterns shown in Fig. |^a) and (b); the pattern re¬ 
peats itself with a period A 2 ; 3.5. The introduction 

of gain/loss to the waveguides gives rise to a chiral local 
structure; two snapshots of such local patterns are shown 
in Fig. I^c) and (d) for a T^T-symmetric kagome lattice 
with 7 = 0.78. The dynamical evolution of the local 
structures is more complex than that of a passive array. 
Again, while maintaining a Cs symmetry, the chiral pat¬ 
tern rotates about the beam axis. This is illustrated in 
Fig.j^b) which shows the angle ©3 as a function of prop¬ 
agation distance z for a narrow beam excitation. The chi¬ 
ral nature of the local structure also manifests itself in 
the beam dynamics: the moving average of 63 deviates 
significantly from zero especially at large 2 :. 

To summarize, we have presented a T^T-symmetric 
kagome photonic lattice. The complex waveguide array 
supports an exact T^T-symmetric phase for gain/loss pa¬ 
rameter below a finite threshold. The eigen-spectrum in 
this phase contains two nearly flat bands inherent from 
the underlying kagome structure. The proposed complex 
lattice possessing a full Cs symmetry can be derived from 
placing T^T-symmetric dimers at either the kagome sites 


or the honeycomb edges. The proposed structure proba¬ 
bly is the simplest complex lattice that respects both the 
VT and C 3 symmetries. We have also uncovered oscilla¬ 
tory power rotation and long-lived chiral local structures 
in the beam dynamics of the complex kagome lattice. 
The presence of nearly flat bands in a T^T-symmetric 
phase could lead to other novel phenomena when dis¬ 
order and/or nonlinearity are taken into account, which 
will be the topics of future investigation. 
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